New CMB Power Spectrum 
Constraints from MSAMl 
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mic microwave background (CMB) at angular scales 
nea r 0?5. The first tw o flights of MSA Ml, reported 
in (pheng et al. 1994D (MSAM92) and ( |Cheng et al 



G. W. Wilson^-^. L. Knox^ S. Dodclson^. K. Goblet 1996| ) (MSAM94), observed overlapping fields on the 
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sky and demonstrated the repeatability of the mea- 
surement. A detailed comparison, showing consis- 
tency between these two fiights, was reported in (In- 
man et al. 1997]) and (|Knox et al. 199"8| ). A 



third flight (Cheng et al. 1997) (MSAM95) mea 



Abstract 

We present new cosmic microwave background (CMB) 
anisotropy results from the combined analysis of the 
three flights of the first Medium Scale Anisotropy 
Measurement (MSAMl). This balloon-borne bolo- 
metric instrument measured about 10 square degrees 
of sky at half-degree resolution in 4 frequency bands 
from 5.2 cm~^ to 20 cm~^ with a high signal-to- noise 
ratio. Here we present an overview of our analysis 
methods, compare the results from the three fiights, 
derive new constraints on the CMB power spectrum 
from the combined data and reduce the data to total 
power Wiener-filtered maps of the CMB. A key fea- 
ture of this new analysis is a determination of the am- 
plitude of CMB fluctuations a,t i ^ 400. The analysis 
technique is described in a companion paper (Knox 



1999) 



Subject headings: balloons — cosmic microwave back- 
ground — cosmology: observations — infrared: ISM: con- 
tinuum 

1. Introduction 

The Medium Scale Anisotropy Measurement (MSAM) 
is a balloon-borne telescope and bolometric radiome- 
ter designed to measure the anisotropy in the cos- 
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sured a nearby region of sky using the same ob- 
serving method. This increased the experimental 
sky coverage and sensitivity to the CMB anisotropy 
power spectrum. A second version of this instrument 
(MSAM II) with complementary frequency coverage 
has since been flown. This data set is still being ana- 
lyzed. 

2. Instrument and Observations 



The MSAMl instrument is described in (Fixsen 
et al. 1996). We give a summary here. The ac- 



tively pointed gondola is composed of a 1.4 m off- 
axis Cassegrain telescope with a multimode bolo- 
metric radiometer. A three-position chopping sec- 
ondary throws the frequency independent '--^0?5 pri- 
mary beam ±0?7 tangent to the local horizon at 
2 Hz. The four spectral channels at 5.7, 9.3, 16.5, 
and 22.5 cm"-'^, each have bandwidth of ^ 1.5 cm~^. 
The detectors' outputs are synchronously sampled at 
32 Hz: 4 times for each of 4 positions of the sec- 
ondary mirror, for a total of 16 samples per chop- 
per cycle. Telescope pointing is controlled with a 
star camera and gyroscope. The configuration of the 
gondola superstructure was changed between the 92 
and 94 flights to reduce possible reflection of ground 
radiation. The improved conflguration remained for 
MSAM95. 

All three flights were launched from the National 
Scientific Balloon Facility in Palestine, Texas. The 



observing method, also described in ( Fixsen et al. 
1999) , is a slow azimuth scan of a region crossing the 



meridian above the North Celestial Pole. For a period 
of 20 minutes the scan center tracks a fixed spot on the 
sky as the earth rotates. Afterwards, an overlapping 
region is scanned. MSAM92 and MSAM94 observed 
at declination 6 = -F82?0, MSAM95 at 6 = -F80?5. 
The fiights observed between right ascensions 14^2 
and 19^.^5. The lower declination of the MSAM95 re- 
quired a faster scan rate because of the increased sky 
motion. The sky coverage of all the MSAMl flights 



is shown in Fig. n^. 




Fig. 1. — Locations in the data file (flat) sky coordinates 
for the 1010 MSAMl points. The boxes show the relative 
twist of the beam-pattern during the observation. The 
upper row of points come from the overlapping flights, 
two years apart, of MSAM92 and MSAM94. The lower 
row are the MSAM95 points. 



3. Data Reduction 

The data from each of the three flights of MSAMl 
are independently reduced in the same manner. A 
detailed discussion of the analysis for each of the three 
flights is available in (pheng et al. 1994, Cheng et al. 



1996, Cheng et al. 1997) respectively. We outline the 



process here. 

1) Spikes caused by cosmic rays are removed from 
the time stream by a filtering and peak detecting tech- 
nique which results in the deletion of 5% to 10% of 
the data. Samples are also lost due to spurious electri- 
cal pickup and telemetry dropouts. For each of these 
cuts, a full chopper cycle is deleted. The total loss is 
between 10 and 30% of the raw data. 

2) The detector time streams are demodulated in 
two ways - each resulting in an independent instru- 
mental beam pattern and corresponding instrumen- 
tal window function. If Tl, Tc, and T/j are the 
sky temperature at the left, center, and right posi- 
tion of the beam during a chopper cycle, the single 
difference demodulation is T/j — T^, making an an- 
tisymmetric beam pattern, while the double differ- 
ence is Tc — {Tl + Tfi)/2, making a symmetric beam- 
pattern. Optimum weighting for the demodulations 
are determined from Jupiter observations. The in- 
strument noise is uncorrelated between between the 
two demodulations. 

3) The data are calibrated using scan and raster 



observations of Jupiter. The brightness temperature 
of Jupiter is reported in (Goldin et al. 1997). Of 
the two models presented in that paper, we use the 



temperatures based on the "Rudy" model ( Rudy et al. 



1987). The error in the calibration is estimated to 



be 5%, dominated by the uncertainty in the Jupiter 
temperature. 

4) The Jupiter raster observations, performed dur- 
ing each flight, are the basis for the high fidelity de- 
termination of the beam pattern for each demodula- 
tion. Beam-pattern uncertainties are dominated by 
arcminute-scale pointing uncertainties. 

5) The estimate of the instrument noise is deter- 
mined from the variance in 100 s segments of the de- 
modulated data after the removal of a slowly drifting 
offset. The offset ranges from 1 to 6 mK in MSAM94 
and MSAM95 with an offset of 10 mK R-J in all chan- 
nels in MSAM92. The drift in the offset is smaU com- 
pared to its value. Because the removal of the offset 
correlates the noise on time scales longer than the 
detector time constant, the remainder of the data re- 
duction incorporates the full noise covariance matrix. 

6) The data are binned according to both the 
position on the sky and the twist of the demodu- 
lated beam pattern during each complete chopper cy- 
cle. The bin size for the twist dimension is deter- 
mined by defining a "Binning Degradation Factor," 
BDF = ^J &■'■ + (S^) jcj where a is the estimated in- 
strument noise and 



dk is an esti- 



(<52)=/(|Bi(^)-B2(^fc)|2Ccdn.(^) 

mate of the expected error in the estimate of the sig- 
nal due to the twist bin size. ((5^) is determined us- 
ing the standard cold dark matter correlation function 
C'cdm(^) convolved with beam-patterns, Bi(k) twisted 
with respect to each other by the rotation matrix A. 
A similar construction is used to define the BDF for 
the spatial binning. The BDF can be thought of as 
the factor by which the sensitivity of the data set is 
decreased due to the choice of bin size. The bin sizes 
are chosen to hold the BDF to values less than 1.1. 
This results in 5° twist bins and 14'bins in sky posi- 
tion. 

7) The calibrated data are analyzed to provide 
measurements of brightness in the four spectral chan- 
nels as a function of bin. The linear combination of 
the spectral channels which minimizes the sensitivity 
of galactic dust foreground and matches the signature 
of a CMB thermal fiuctuation over the spectral range 
of the instrument channels is found and an estimate of 



CMB anisotropy and dust optical depth for each bin 
is produced. This is done by fitting the data for each 
bin in the four channels to a two parameter model of 
sky and dust. 

4. Comparison of MSAM92 and MSAM94 

The overlapping regions of MSAM92 and MSAM94 
(see Fig. P are used to compare their estimated sky 
signals. This can place a limit on how much of the 
signal could be attributed to instrumental artifacts or 
other local contamination. While straightforward in 
principle, a simple comparison is not possible despite 
the large degree of overlap. The beam centers for each 
sample do not line up perfectly and because of the 
twist dimension in the binning, there are few bins that 



are populated in both flights. In (Inman et al. 1997) 
the bin size was expanded over the criterion in the 
previous section and those bins with sufficient data 
were differenced. With rather reduced sky coverage, 



(Inman et al. 1997) found no signal in the differenced 
data. 

An alternative procedure for comparing the two 



measurements has been previously reported in (Knox 
et al. 1998). Here, an assumed power spectrum for 



the, CMB fluctuations is used to make a prediction of 



thel most likely signal in the 1992 data sot, given the 
1992 pointing information. This is compared to the 
most likely value of the signal in the same 1992 data 
set but given the 1994 pointings. This "most likely" 
signal is determined by applying a Wiener fllter to the 
data. See ( Knox et al. 1998 ) for details. In Fig. || we 
see that the two data sets predict very similar signals 
for the 1992 data set using either the 1992 or 1994 
data. 

We quantify "very similar" by use of the likelihood 
ratio statistic. The two hypotheses are 1) the sig- 
nals are correlated as one would expect (given the two 
sampling strategies and an assumed power spectrum) 
and 2) the signals are uncorrelated between data sets. 
We use the natural log of the likelihood ratio statistic, 
which is a quadratic operator on the dat a denoted by 
(3 f\ For the 1992 and 1994 data sets ( |Knox et al 



1998) has found (3 = 12.8, which means that hypoth- 
esis 1 is e^^'* times more likely than hypothesis 2. A 
frequentist interpretation of P is given by calculat- 
ing the expected mean and standard deviation of the 
statistic under the different hypotheses. The result is 



^Also see ( Fegmark 1998) on the optimization of quadratic com- 
parison statistics. 
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Fig. 2. — Most likely signal in 1992 data set, given the 
COBE-normalized standard CDM power spectrum and 
the 1992 data (vertical lines) or the 1994 data (horizon- 
tal lines). The shaded area is the 68% confidence region. 
Single-difference (or "2-beam" ) data in top panel, double- 
difference (or "3-beam") in the bottom panel. 



15.0±4.1 (hypothesis 1) and -58.4 ±27. 4 (hypothesis 
2). This analysis is in agreement with (Inman et al. 



1997) that it is extremely unlikely that the data sets 



are caused by a signal that is uncorrelated between 
experiments. Based on these analyses, we conclude 
that the signal comes from the sky and not from the 
instrument or local environment. 

5. Likelihood Analysis 

The data set from the three flights of MSAM has 
been reduced to 505 measurements of the CMB sky, 
for each of the two demodulations. We model this 
data, d as due to signal and noise 



Si + n. 



(1) 



where i runs from 1 to 505 over the single-difference 
demodulation and from 506 to 1010 over the double- 
difference demodulation, and the signal, s, is related 
to the true temperature field, T, by 



Si = B{x — Xi)T{x)dx. 



(2) 



Here, B{x) is the (single or double difference) beammap, 
and Xi specifies the pointing. We assume that both 
the signal and noise are Gaussian-distributed with 



zero mean with covariance matrices which we denote 
by Sij = (siSj) and iVy = {nirij). 

The noise covariance matrix, N, is block-diagonal 
with each block representing the noise correlations of 
a single demodulation from a single flight. The noise 
covariance matrix is singular due to the independent 
offset removals from each of the three flights (two each 
in MSAM92 and MSAM94). This constraint must be 
explicitly projected out of the data which we do with 
an SVD inversion of N. 

The signal covariance matrix is linearly related to 
the angular power spectrum of the temperature field, 
Ci . The likelihood of this power spectrum, given the 
data, noise matrix and our assumptions of Gaussian- 
ity is 

L{Ci) = — , , (3) 

where d is the n-element vector of observations and 
C = [S{Ci) + N] is the (n x n) covariance matrix 
of the observations. We use this likelihood to place 
limits on the power spectrum of fiuctuations. 

For this analysis we parameterize the theoretical 
signal covariance matrix, S, with the power spectrum, 
Ci = l{l + l)C;/(27r), broken into bands denoted by B 
so that 

Ci^J^Xs^n^B, (4) 

B 

where Cb denotes is a flat power spectrum within 
band B with amplitude Cb- That is, 



^B = -Z ^B- 

III 

The sum runs over the bands in /-space with 



Xe 



1 : U{B)<l<ly{B) 
: otherwise 



(5) 



(6) 



This parameterization of Ci is completely general and 
its usefulness will become apparent below. 

The calculation of the likelihood requires the inver- 
sion of the (n X n) covariance matrix C . It has been 



shown in ( Bond 1994 ^ Tegmark, Taylor, and Heaveng 



19c(6|, |Bunn and White 1994 |Bond and Jaffc 1997| ) 



that a substantial reduction in the rank of C can be 
achieved by working in the signal-to-noise eigenmode 
basis. This is true even in a high signal-to-noise case 
like that of MSAMl. For this data set, we achieve 
a compression by a factor of 1.8 in the rank of C by 
ignoring modes with signal-to-noise ratio of less than 



0.03. Working in the signal-to-noise eigenmode ba- 
sis has the added benefit of automatically projecting 
out the eigenmodes associated with the offset removal. 
Thus, only one initial SVD of the covariance matrix is 
required (to zero the infinite eigenvalues). Inversions 
of the covariance matrix in the signal-to-noise eigen- 
mode basis are then done using faster methods such 
as Cholesky decomposition. 

6. The Flat Band-Pov^rer 

As has been done previously for each of the data 



sets individually (Cheng et al. 1994, Cheng et al. 
1996, Cheng et al. 1997), we calculate "flat band- 
powers" for each demodulation. That is, we as- 
sume the entire power spectrum is flat with amplitude 
Ci = Cb and calculate the likelihood of this amplitude. 
Table gives the flat band-powers (maximum likeli- 
hood values of VCs) for the three flights of MSAMl 
for the single and double difference demodulations. 
The error bars indicate where the likelihood falls to 
e~^" of the maximum. 



Flight 


Single Diff. 


Double Diff. 


MSAM92 


48 ±11 


54 ±10 


MSAM94 


35 ±6 


45 ±9 


MSAM95 


51 ±7 


56 ±7 


all three 


47 ±5 


53 ±5 



Table 1: Flat Band-power Estimates for MSAMl in 
/iK 



7. Radical Compression 

Flat band-powers, together with the diagonal parts 
of the window function matrix, often simply called the 
window function, have traditionally been the main 
results of CMB experiments. When taken together, 
they are the raw ingredients for constraining the 
power spectrum and cosmological parameters. 

The parameters, Op, are found by minimizing the 
X^ where 



X 



^C^lBiCi{ap) 



Cb?Io 



Si 



(7) 



where B runs over different data sets, Cb and ob 
are the band-powers and their standard errors respec- 
tively, and ]bi is a filter which, when summed over 
the power spectrum, Cj, gives the theoretical predic- 
tion for the band-power, Cb- The filter is usually 



constructed from the familiar window function, Wbi 

Wbi/1 



Ibi — 



iY^WBiiiy 



(8) 



We call the filter given by this equation the window 
function filter. The parameters Up could be cosmolog- 
ical parameters (e.g., fi^, ^a, Hq, etc.) or parameters 
from a phenomenological power spectrum. 

7.1. Problems With Flat Band Powers ... 

Using flat band-powers as a form of radical com- 
pression has the following drawbacks: 

1. The actual sky power spectrum is not flat. 

2. The expectation value of the band power is not 
given by summing the window function filter 
over the power spectrum, and thus the window 
function filter should not be used in Eq. 0. The 
expectation value is only given by this sum in 
the limit that the data points have no signal 
correlations. 

3. The method provides no estimate of the corre- 
lation between the errors in the estimates of Cb 
from different demodulations. 

4. The constraints on the parameters are not Gaus- 
sian, even though this assumption is implicit in 
the x^ minimization. 

Problems H and H are well known deficiencies of 
the bandpow er approach. Problem ^ has b een em- 
phasized in ( [Bond, Jaffe, and Knox 1998b ), where 
an approximate solution was given. Here we focus 
on problem El which has been discussed previously in 
( Bond, Jaffe, and Knox 1998a ). We illustrate the po- 
tential severity of the problem with an extreme exam- 
ple. Consider a total power mapping experiment with 
angular resolution of FWHM = 30' which has mea- 
sured a 5° X 5° patch of the sky. The window function 
filter for this experiment is /; oc Wi/l — exp{—Pa'^)/l 
where tTf, ~ FWHM/2.355. Note that this filter peaks 
at I = indicating that the experiment is most sensi- 
tive to flucutations on very large angular scales. How- 
ever, the data set is not actually sensitive to the lowest 
spatial frequencies at all. The problem lies in having 
ignored the off-diagonal terms. The filter function ac- 
tually makes sense only if the data points are all far 
apart on the sky so that Sij is diagonal. For the exam- 
ple given, correlations between the points on the sky 



are making the data set insensitive to fluctuations on 
large scales. Because using the diagonal component of 
the window function to define the filter ignores these 
correlations, we get a nonsensical result. For most 
actual data sets, the problem is not quite so severe 
but this example illustrates the potential pitfall. 



7.2. 



and Solutions 



Solutions have been found to all of these problems 



(Bond, Jaffe, and Knox 1998a, Bond, Jaffe, and Knox 



1998b, Knox 1999). Here we briefly review them. 



Problem n^ can be solved by breaking the power 
spectrum into several bands as in Eq. 0, and then 
finding the amplitudes of these bands, Cb, that max- 
imize the likelihood. We find this maximum by it- 
erative application of a quadratic estimator, as has 



been done f or COBE/DMR Qcnnctt et al. 1996D and 
Saskatoon ([Netterficld et al. 1997 ) data in (Bond, 
Jaffe, and Knox 1998a|) and on simulated MAP data 



in ( Oh, Spergel, and Hinshaw 1999| ). By calculating 
the covariance matrix of the set of Cb we also solve 
problem y. 

Because physical power spectra are not actually 
fiat across these bands, we need a means of taking 
a general power spectrum, Ci, and turning it into a 
prediction for Cb- In other words, we need to be able 
to calculate the expectation value of Cb, (Cb), under 
the assumption that the power spectrum is Ci. This 
relationship is specified by the filter function. 



{Cb) — IbiCi 



(9) 



Taking Eq. ^ as the definition of the filter, ( pK^nox| 
1999| ) has shown how it can be calculated from the 
signal and noise covariance matrices and the deriva- 
tives of 5' with respect to C;. Taking into account 
all off-diagonal terms, this prescription for fsi solves 
problem 0. To distinguish it from the usual practice 
of simply using the window function filter, we call this 
the minimum-variance filter. They are identical only 
in the limit of no signal correlations. 

We could remove the need for filters by making the 
bands very narrow since sufficiently narrow bands en- 
sure that the sensitivity to each Ci within the band is 
approximately independent of Z. However, making the 
bands too narrow would increase the non-Gaussianity 
- exasperating problem H - because the likelihood of 
more tightly constrained broad bands is better ap- 
proximated by a Gaussian. Therefore, the bands must 
be broad enough to have significant constraints on 



their amplitudes. For MSAM, this condition makes 
the bands sufficiently broad that the sensitivity to Ci 
varies significantly across the bands, necessitating the 
use of a separate filter for each band. 



Finally, if we adopt the (Bond, Jaffc, and Knox 
19!! 8a) prescription for problem W, which requires cal- 
culation of a "log-normal offset" , xb, for each Cb, we 
have solutions to all four problems. Although these 
solutions are not exact, they do represent a significant 
improvement over the usual flat band-power method. 

It is not necessary to break the power spectrum 
into bands to obtain parameter estimates from the 
observations. However, this approach aids the com- 
parison of different experiments with a minimum of 
theoretical assumptions, as well as easing the compar- 
ison of experimental results with theory. By following 
the above procedure for power spectrum estimation, 
the full weight of an experiment is made available in 
an easily-tractable form for the kind of parameter es- 
timation outlined in Eq. [^ 

7.3. The Application to MSAM I Data 

The MSAMl data sets are a prime example of 
the limitations of the flat band-power method. The 
MSAMl data have high signal-to-noise and are heav- 
ily sample-variance limited when using standard esti- 
mators of the flat band-power. We now use the radical 
compression methods outlined above to probe regions 
of ^-space ignored by our previous reduction to flat 
band-powers. 

The difference between the minimum-variance fil- 
ters, /s;, and the window function filters, is shown in 
Fig. g. The plot shows the filters for each demodula- 
tion of the 3 year data for a single band covering all I. 
Notice that the minimum-variance filters show more 
response at high I than the window function filters. 
This is due to the fact that the dense sampling and 
high signal-to-noise ratio of the data set yield infor- 
mation on angular scales smaller than the beam size. 
Again, this information is in the off- diagonal com- 
ponents of the covariance matrix - underscoring the 
need for experiments to track the full noise covariance 
matrix in the data reduction. 

We plot the minimum-variance filters for the in- 
dividual demodulations only to make the point that 
they are not equal to the window-function filter as has 
often been assumed in the past. In the analysis we 
describe below, we do not treat the double-difference 
data sets and single-difference data sets separately; 




f. 



Fig. 3. — The minimum-variance filters (solid lines) and 
window function filters (dashed lines) for the single and 
double difference demodulations. 



the very significant correlations between them are in- 
cluded. 

For this analysis, we break up the ^-space coverage 
into three wide bands and allow Cb to vary in each. 
In line with the discussion above, we choose the three 
bands such that each has enough weight to produce 
an interesting constraint on the power spectrum. The 
Z-ranges for the three bands chosen are 39-130, 131- 
283, and 284-806. 

In Fig. we plot the 3 power spectrum estimates, 
Cb from the combined three years of data. The cen- 
tral points are located in Z-space at l^s which we de- 
fine as the average value of I over the filter function 
for that band. We include a horizontal bar from Z_ to 
Z+ . For the middle band, these are taken to be simply 
the beginning and end of the band (L = Z< = 131, 
Z+ = Z> = 283). For the far ends of the two outer 
bands, we take them to be where the filter falls to 
e~^''^ of maximum. Similar results (dashed lines) are 
achieved by analyzing a total-power map of the CMB 
temperature which is constructed from the demod- 
ulated data. We take the good agreement between 
these two nearly independent analysis techniques as 
strong proof that we have calculated the complicated 
signal covariances correctly. We will discuss this pro- 
cedure in section ^. 

While the Cb are not independent, their correlation 
coefficients are fairly small. The correlation between 
bands 1 and 2 is —0.18, between bands 1 and 3 is 
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Fig. 4.— The MSAM band-power estimates. The 
sohd hnes give the estimates of the power in the three 
bands calculated directly from the demodulated data, 
and recorded in Table |[ Similar results (dashed lines) 
are achieved by analyzing a total-power map of the 
CMB temperature which is constructed from the de- 
modulated data (see section 0). 



-0.024 and between bands 2 and 3 is -0.29. The 
error bars shown in Fig. H are the result of marginal- 
izing over the power in the other bands. Under the 
assumption that the other bands are fixed, the error 
on the band in question is less than 5 % smaller. 

The power spectrum estimates, Cb, their weight 
matrix, Fui , filter functions, fsi, as well as log-normal 
offsets, xb (see (Knox 1999)) are available at 
http://topweb.gsfc.nasa.gov and also at 



http://www.cita. utoronto.ca^knox/radical.htmj which 
includes similar information from other CMB anisotropy 
data sets. 

8. CMB Maps 

A useful check of our power spectrum results can be 
made by analyzing a map made from the demodulated 
data as opposed to directly from the demodulated 
data as we have done above. We begin constructing 
this map by recognizing that Eq. 1 and 2 can be 
combined and rewritten in matrix form as 



d = BT 



(10) 



With the assumption that the noise is Gaussian, with 
covariance matrix, N ^ the most likely value of T, given 
the data, d, is that which minimizes the x^: 



X 



{d- BT)N-^{d~ BT). 



(11) 



This minimum, which we denote by T, is given by 

f = NBN-^d. (12) 

This estimate of T will be distributed around the true 
value due to noise, where the noise covariance matrix 

is 



TV EE< (T - T){T - T) >= {B'^N-^B 



(13) 



Table 2: Power Spectrum Estimates from MSAMl 



l- 


IcS 


1+ 


^Cb (/iK) 


39 


84 


130 


351- 


131 


201 


283 


49±^" 


284 


407 


453 


471^ 



This map can be analyzed in the same manner as 
the demodulated data, with the advantage that the 
signal covariance matrix is now very simple to com- 
pute. Previously, calculating the signal covariance 
matrix required doing a four-dimensional integral for 
every covariance element. In this new "map basis," 
the signal covariance matrix simplifies to 



v^ 21 

I 



1 



An 



F,(cos(%))Ci. 



(14) 



The price to pay for this simplicity is that the noise 
covariance, iV, is very complicated. We have done this 



analysis as a check of the calculations in section [t!^ . 
The results are shown in Fig. |[ The agreement is a 
strong argument that we have made no errors in what 
is a fairly elaborate and difficult calculation. 

The map, T, is extremely noisy and not visually 
useful. We can greatly reduce the noise by Wiener fil- 
tering, (e.g.jBunn, Hoffman, and Silk 1995, Tcgmark 
et \il. 1997|, |Knox et al. 199^ ). The~Wiener fiher 



produces the most likely T, given not only the data, 
but also an assumed power spectrum for the signal. 
The Wiener filtered maps are shown in Fig. a. 
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Fig. 5. — A map of the three years of data. Top region 
was covered by the 92/94 flight; bottom by 95 flight. 



9. Discussion 

Because the third band is derived by making mea- 
surements on scales on order of the beam size, we 
must ask what sensitivity the amplitude of Cb has to 
the beam shape. For example, if the band sensitivity 
results from high frequency fluctuations in the mea- 
surements of the various Bi(k), the estimation of the 
amplitude of Cb would be suspect to the errors in de- 
termining the beammap. We address this question by 



performing a number of analyses of ihe three years of 



datfe The first analysis (leadmg lo the quoted val- 
ues of Cb here) is done using the six beammaps mea- 
sured in the three flights. That is, MSAM92 data 
goes with the MSAM92 beammaps, MSAM94 data 
with MSAM94 beammaps, and MSAM95 data with 
MSAM95 beammaps. The second analysis is done 
using the beammaps measured during the MSAM92 
flight for all three years of data. This is followed 
by repeating the analysis with the beammaps from 



MSAM95 0. We find that after accounting for the 
normalizations of the different data sets, there is no 
evidence for the third band being sensitive to the 
beammap choice. 

Reanalyzing the entire data set using beammaps 
measured in flight from raster observations of Jupiter 
for each of the three different years the experiment 
flew, is taken to be the most pessimistic estimator of 
the effect of the beam on the third band. The differ- 
ences between the beammaps include all the statisti- 
cal errors of the beammaps, any errors in the raster 
observations themselves, and any changes introduced 
by the complete rebuild, realignment, and refocusing 
of the optical system and instrument configuration. 

To place our estimates of the power spectrum in 
context, we plot them with the predictions of several 
theoretical models as well as a fit of the power in 11 



bands to all available data (based on (Bond, Jaffe 



and Knox 1998b| ) including the previously published 



MSAMl points. 
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Fig. 6. — The dark error bars are the power spectrum 
constraints from the 3 year MS AM data set. The 
light error bars are the result of a fit of the power 



in 11 bands to all ava ilable data (based on (Bond 
Jaffe, and Knox 1998b| ) including the previously pub- 



lished MSAMl points. The curves are standard CDM 
(solid), a ffat Lambda model (dotted), and an open 
model with ficurvaturc — 0.6 (dashed). 



^The MSAM92 and MSAM94 beammaps are similar enough to 
be swapped with no change in estimated signal. 



10. Summary 

We have calculated new power spectrum estimates 
from the combined three flights of the MSAMl in- 
strument. The analysis technique used is an improve- 
ment over the standard flat band-power aproach and 
includes all correlations in the data. In addition to 
power-spectrum estimates and their error covariance 
matrices we have also provided the log-normal off- 
sets and minimum- variance filters in order to improve 
"radical compression." The analysis yields a strong 
constraint on the power spectrum at I ~ 400, broad- 
ening the /-space coverage of the experiment into a 
theoretically very interesting region. 
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